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1. Introduction 

The study of processes involving the production of high transverse momentum [pt) QCD 
jets will form an integral part of the LHC physics programme. In this light it is clear that 
the understanding of properties of jets will impact the accuracy of studies concerning new 
physics as well as traditional QCD studies concerning for instance the extraction of parton 
distribution functions. 

While the subject of jet physics is far from new and has seen steady progress since 
the pioneering work of Sterman and Weinberg Q, the LHC provides a challenge of un- 
precedented complexity owing to the initial state hadronic environment and the very high 
beam energy. The QCD tools at our disposal have therefore to be developed and adapted 
to meet this challenge. 



- 1 - 



The most important tool in the context of accurate jet studies at the LHC will remain 
QCD perturbation theory. Reliable perturbative estimates will constitute the only first 
principles (model-independent) information available from QCD theory and hence accurate 
perturbative results for various processes of interest are imperative. The advent of new 
infrared and collinear (IRC)-safe jet algorithms such as SISCone Q and the anti-A;^ Q 
jet algorithms to complement existing IRC-safe algorithms such as the kf [Q, ^ and the 
Cambridge/ Aachen algorithms Q is therefore an encouraging and crucial development. 

Perturbative aspects apart, another serious stumbling block for jet studies at the LHC 
is the role of non-perturbative effects. Here one is dealing not just with the hadronisation 
of partons but also the underlying event arising from beam remnant interactions, multiple 
hard interactions and the issue of pile-up. Thus while in an ideal world one may hope to 
construct clean signals for new physics, for instance mass peaks signifying the production 
of new particles, in the real hadron collider environment such peaks can be significantly 
smeared by non-perturbative QCD effects such as hadronisation and the underlying event, 
in addition to perturbative radiation. A better understanding of the size and role of non- 
perturbative effects as a function of the parameters used in jet studies (such as jet size) is 
hence imperative. 

While neither hadronisation nor the underlying event can be calculated within pertur- 
bation theory one could argue that one may have more control on the former effect. The 
reason for this argument is that the hadronisation of partons that constitute a high-pj jet 
is a phenomenon that has been met and handled before, for instance at LEP, HERA and 
the Tevatron. There has therefore been considerable development of hadronisation models 
in Monte Carlo event generators such as PYTHIA Q and HERWIG 

Also significantly, in addition to Monte Carlo simulations, analytical techniques based 
on renormalons have yielded important information on hadronisation corrections in the 
context of event-shape variables at LEP and HERA . The pioneering phenomenological 
studies in this regard were carried out by Dokshitzer and Webber Using a model for a 
universal infrared- finite coupling they were able to compute the 1/Q power corrections to 
event-shape variables at LEP. The ensuing predictions involved just a single universal non- 
perturbative parameter, the moment of the coupling over the infrared region. Extracting 
the value of this moment from data on one variable it thus became possible to predict 
the 1/Q corrections to a large class of event-shape variables. Similar success was also met 
in the description of Breit-frame current-hemisphere event-shape variables at HERA |12| 
in spite of the presence of an initial-state proton that may potentially have impacted the 
power corrections. 

The study of hadronisation corrections at hadron colliders, based on analytical tech- 
niques, is however still in its infancy. The main deterrents to such studies have been for 
instance the more complicated colour topology of multijet hard processes at hadron col- 
liders and the competing presence of the underlying event. The issue of the role of colour 
in influencing the hadronisation corrections was explained in ref. |13| in the context of 
away-from-jet energy flow. 

As far as the role of the underlying event compared to hadronisation is concerned 
an interesting result was derived in ref. |14]. There it was observed that if one studies 
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the hadronisation contribution to the jet energy (specifically the change in pt, dpt of a 
jet due to hadronisation) using the aforementioned analytical techniques, one obtains a 
singular dependence^ on the jet radius R with the hadronisation varying as —1/R. The 
underlying event on the other hand makes a contribution to the jet pt varying as (see 



also refs. |16, The possibility thus emerges of a parametric separation of the two 



components of non-perturbative QCD with hadronisation being the dominant effect for 
jets at small R. Moreover the leading small- i? behaviour was seen to arise due to coUinear 
soft gluon radiation from the triggered hard parton jet. The colour factor associated to 
collinear radiation is merely the colour charge of the emitting hard parton, which means 
that in this limit the complex colour topology of the whole event is essentially irrelevant. 
Knowing the i?-dependence of perturbative and non-perturbative corrections then allows 
for the possibility of deducing optimal R values adapted to different studies at hadron 
cohiders 0. 

Perhaps most interestingly a link was also made in ref. |Q] between the magnitude 
of the 1/R hadronisation correction to the jet pt and that of 1/Q power corrections to 
event shapes such as the thrust in e~^e~ annihilation |11]. The coefficient of the 1/R 
contribution is a perturbatively calculable number times a non-perturbative factor which 
is the same coupling moment as enters the event shape predictions. Carrying out the 
relevant perturbative calculation at the single-gluon level the coefficient of the 1/R 
correction to the jet pt for a quark jet is seen to be one-half of the coefficient of the 1/Q 
power correction for the e~^e~ thrust variable as computed in ref. []TT| . 

The predictions of refs. |jll| and iQ are based on one-loop 0{as) perturbative calcu- 
lations. At this level the emission of a soft gluon with transverse momentum kt ~ Aqcd 
is associated to non-perturbative hadronisation contributions to the transverse momentum 
Pt of a hard jet or the event shape under consideration. Physically this is due to the fact 
that the running coupling associated to such an emission is as{kt) which is perturbatively 
undefined at kt ~ Aqcd- Replacing the unphysical perturbative coupling in this region by 
a finite and universal non-perturbative extension^ leads to the predictions for hadronisation 
corrections reported in refs. [11| and [p!^]. 

In the context of event shapes however it was shown by Nason and Seymour that one- 
loop calculations involving single gluon emission are inadequate to predict the coefficient 
of the 1/Q power corrections |19|. The argument of the running coupling which we took 
to be kt can in fact only be reconstructed at the two-loop level where one has to account 
for the decay of the emitted gluon into offspring gluons or a quark/anti-quark pair. For 
variables such as event shapes and some observables involving jet definition such as jet pt, 
one is not free to inclusively integrate over the decay products since the observable depends 
on the precise details of the decay kinematics. Thus in these cases the replacement of a 
gluon decay with a parent virtual (massive) gluon, with a running coupling dependent on 
the parent virtuality or kt, is not actually correct at the level of power corrections. One 
must then address the gluon decay and identify the correction needed to the single gluon 



similar result was obtained by Korchemsky and Sterman for cone jets in e^e annihilation (see ref. 

iii)- 

^The treatment of the running coupling was formalised via the dispersive approach adopted in ref. [hsl. 
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results quoted for example in refs. [11| and |14]. 

For the case of the thrust and other commonly studied e^e~ event-shape variables cal- 
culations accounting for gluon decay have already been performed some time ago |2^, ^ 
and yielded a simple result: taking into account gluon decay resulted in a universal (event- 
shape-independent) factor multiplying the single massive gluon results - the "Milan" fac- 
tor. This was also seen to be the case for DIS event shapes |22| and further confirmed by 
the results of ref. p^]. One of the key ingredients in obtaining a universal Milan factor 
was the fact that all the event-shape variables considered are linear in transverse momenta 
of soft emitted particles |21]. The jet pt variable is however not so simple. For almost all 
the commonly used jet algorithms it turns out that the 5pt due to gluon emission is not 
simply a linear sum over the contributions of individual gluons, and that the contribution 
of a given gluon depends on whether or not it is clustered or combined with other emis- 
sions. This effect is absent for the anti-fej algorithm where soft gluons essentially cluster 
independently to the hard emitting parton and the result is the universal Milan factor 

In the present paper we establish that for jets defined in, for instance, the kt algorithm 
(which we treat explicitly) the resulting rescaling factor is not the universal Milan factor. 
However we prove that it is a calculable factor of the same order, which we compute for the 
kt algorithm. Thus it is now indeed possible to link the results for hadronisation corrections 
to jet Pt with those obtained for event shapes after including the Milan factor. As a result 
of our calculations the relation quoted in ref. |14| between the jet pt hadronisation and the 
thrust power correction changes for all but the anti-fet algorithm. The hadronisation cor- 
rections to jet Pt become algorithm-dependent rather than truly universal. They however 
remain universal in the more crucial sense that the result is always a perturbatively calcu- 
lable number multiplying a universal non-perturbative coefficient. The only complication 
is that the perturbative calculations have to be at the two-loop level. 

The current paper is organised as follows. In section ^ we deal with the impact of gluon 
emission and decay on the jet pt considering both the kinematical effect on the observable 
as well as the dynamics of gluon branching. The naive one-gluon term calculated in ref. |jl^ 
is identified and seen to be accompanied by two distinct correction terms, which following 
refs. |2C, 21] we refer to as the inclusive and non-inclusive corrections. Section ^ deals with 
extracting the \/R behaviour for the naive result (where we recover the answer quoted in 
ref. |l^]) and computing the inclusive correction. Section |^ is devoted to the treatment of 
the non-inclusive term for the kt algorithm and the extraction of the leading 1/R behaviour 
in this term. In section ^ we combine the results to derive an overall rescaling factor to the 
result of ref. |14| for the kt algorithm. We also recover the expected result for the anti-/ct 
algorithm which agrees with the Milan factor found for event shapes |20, 21|. Finally we 
draw some conclusions and discuss prospects for further work. 

We should also remark that in the following sections we make rather heavy use of 



notation and terminology as well as many results common to refs. pO|, 21, WM, since most 



aspects of the calculation here carry over unaltered from those papers. Thus a full technical 
understanding of some ideas used in this article may require a reading of those references. 
The main ideas and the essential calculations should however be clear on a reading of just 
the current paper. 
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2. The jet 6pt and gluon splitting 



Here we repeat the arguments, first detailed in ref. [pc| ], that led to the emergence of 
the universal Milan factor, but do so in the context of the jet 6pt. We start with some 
introductory remarks that explain the simplifications we make which lead us to the final 
results. 



2.1 Preliminary remarks 



In ref. |14] it was noted that the dominant piece of the non-perturbative hadronisation 
contribution to the jet pt at small jet radius R scaled as 1/R. Additional terms 0{R) 
and beyond were also obtained but these had a negligible impact at small^ R. Due to its 
significance we shall concentrate in this article on the dominant 1/R piece of the result. 
It was also made clear in ref. [14| that the 1/R term arose from gluon emission at the 



boundary of the jet, which in the small- i? limit is in the region collinear to the hard parton 
initiating the jet. To be precise the singular dependence on i? is a direct manifestation of 
the collinear singularity present in the splitting of massless partons due to gluon emission. 
The same collinear singularity was also shown to be associated with InR terms in the 
perturbative regime, which we do not discuss further here. 

Given its origin in collinear radiation it is simple to deduce that the 1/R term does 
not depend on the details of the process of which the triggered hard parton (jet) is a part. 
This fact considerably simplifies the calculations involved and rather than considering all 
possible hard emitting dipoles in a given hard process (such as dijet production in hadron- 
hadron collisions treated in refs. |]l^, |l^) one need only consider collinear branching of the 
massless parton corresponding to the triggered hard jet. 

In spite of this obvious generality, to illustrate the calculation and clarify the connection 
with the existing treatment |14| it may be useful to take a specific example. The example 



chosen in ref. |14] was the change in pt due to hadronisation, 6pt, of a jet produced in 
the threshold limit of hadronic dijet production where threshold was defined by the limit 
Pt Then the change in pt from its threshold value due to soft gluon emission 

was computed and the gluon emission was averaged over to obtain the mean value of 



6pt. In what follows we shall repeat the calculation of ref. |14], this time in the collinear 
approximation sufficient to generate the 1/R behaviour but with account of the decay of 
the emitted gluon. 

2.2 Gluon emission kinematics and dynamics 

We shall also be making use of the formulae and notation used in the Milan factor compu- 



tations and first presented in ref. |2C]. To this end let us introduce Sudakov vectors along 



and opposite to the triggered jet direction, taken itself to be at ninety degrees with respect 



^As a matter of fact the 1/R behaviour was shown to hold to a good approximation over virtually the 
entire range of R values up to /? = 1, which can in part be attributed to the smaller coefficients that 



emerged for the finite R pieces in the analysis of 
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to the beam as required by the threshold Hmit: 



P = ^(1,1,0,0), 

P= ^(1,-1,0,0). (2.1) 



We parameterise the emitted parton four-momenta as in ref. |20|: 



ki = aiP + l3,P + ku, (2.2) 

with kt the transverse momentum vector with respect to the jet direction. For massless 
partons we have that fc? = or aifii = kf^/ s. Also with neglect of recoil of the hard partons 
against the soft parton momenta we can identify P and P with the momenta of the final 
state hard partons themselves (recall that we are considering threshold where the hard 
partons are produced back-to-back and at ninety degrees to the beam). 

We shall also need the value of 5pt due to soft parton emissions. Here one can make 



use of the results derived in ref. |14], where it was shown that emitted partons that are 



recombined with the triggered hard parton contribute to the 5pt by lending a mass MJ to 
the jet. Soft partons not recombined with the jet lend a mass to the "recoil jet" involving 
the recoiling hard parton which is not triggered. Thus one has: [14|: 



where the above equation is correct to first order in soft emitted parton momenta or 
equivalently the small jet masses. In terms of the Sudakov variables defined above one 
obtains: 

^1 1 




2J^ 2J^ \ ^ 2 

where the sum runs over all emissions recombined with the triggered hard parton, to form 
the triggered jet, and we neglected all terms quadratic in the transverse momenta of soft 
emitted partons including those that give a tiny recoil to the hard parton momentum p, 
which we thus identified with the Sudakov vector P. Likewise the mass of the recoil system 
is: 

M 2 _ ^ 



2^ 

where the sum is over all emissions not recombined into jet j. We shall use the above 
results in the following subsection. 

2.3 Calculation of (Spt) with gluon decay 

Here we shall consider the change in pt derived above, together with the squared matrix el- 
ement for gluon production and decay, to obtain the average 5pt- We consider the situation 
up to the two-loop O (q^) level and hence we need to account for single gluon emission, 
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virtual corrections to it as well as correlated two parton emission. To be precise we can 
write H^, H: 



('^P<) = V / ^^{'^siO)+47rxik^)}5ptik)+4CF J (^^Y dr2^6pt{h,k2), (2.6) 

where the first term on the right-hand side of the above represents the contribution from real 
massless gluon emission and the virtual correction to single emission denoted by x{kf). The 
second term on the right-hand side is the contribution from correlated two-parton emission 
with the squared matrix element for gluon decay containing both gluon decay into a 
pair of gluons as well as the Abelian contribution from gluon branching to a qq pair. The 
measure dT2 is the decay phase-space and 6pt{ki,k2) is merely the change in pt induced 
by two-parton emission, a special case of the general situation discussed above. 

We should also mention the fact that the above result is in fact a "dressed" two-loop 
result where an infinite set of higher order perturbative graphs are implicitly embedded in 
the gluon propagator. These graphs are, at least in the Abelian limit, essentially self-energy 
or bubble insertions (renormalon graphs - see |^] for a review) that drive the scale of the 
coupling to the gluon virtuality A;^, hence triggering non-perturbative contributions from 
the infrared region k'^ ~ Aqqj-,. Since for a real gluon the virtuality /c^ is zero, one has the 
ill-defined quantity as{0) appearing above (which will eventually disappear), while in the 
gluon decay case the argument of is k'^ = m?. 



Note also that in writing eq. (2.6) we have specialised to the collinear limit by ignoring 
the fact that the soft gluon k can be emitted by several dipoles that form the hard particle 
ensemble, since in the collinear limit the emission is essentially off the triggered hard parton 
alone. We have taken this parton to be a quark and hence the appearance above of the 
colour factor Cp, but for a gluon jet this can straightforwardly be replaced by Ca- 



We now follow all the arguments of refs. |20, We first attempt to treat 5pt as 
an inclusive quantity which corresponds to the treatment in ref. |14]. We then identify 
that there is a non-vanishing correction needed to the inclusive treatment which we shall 
compute. To this end we write: 



Spt{ki,k2) = dpt{ki, k2) - 6pt(ki + k2) + 5pt{ki + ^2), (2.7) 

where the term 5pt{ki + k2) corresponds to an inclusive treatment where the true value 
of 5pt{ki,k2) is replaced by the value of 5pt that would be obtained by considering just 
the emission of the massive parent gluon ki + k2, with {ki + ^2)^ = inn?. The remainder 
Spt{ki, k2) — Spt{ki + ^2) leads to a correction term that we evaluate in due course^. 



*It should be easy to see that this correction term is in general non-zero. If one considers for example 
a gluon decay where only one of the decay products say ki gets recombined with the jet, then one has 
8pt{k\, ^2) = — (Qf2 + Pi). Assuming that the parent gluon is outside the jet the inclusive contribution 
is Sptiki + k2) — — v's/2 Q — — {ai + a^) which differs from 5pt{k\ , fe). 



-7- 



Substituting (^j) into (^^) we can write: 

IT J nkf a 



+ '^CfJ [^) dr2^6pt{h + k2) + {dptr , (2.8) 

where we separated the non-inclusive correction (denoted ni) from the inclusive approxi- 
mation such that: 

{6ptr = 4CF I (^^ydr2^{5pt{h,k2)-6pt{h + k2)). (2.9) 

For the subsequent discussion we will need the results below for the two-parton (gluon 
decay) phase-space |2^, 21 1: 

fi, I, \ ^daiCpkti da(Pkt(fq 
dr2{ki,k2) = = ^(1 - z)dz , (2.10) 

a, TT a TT TT 

i=l 

with a = oi + a2 the Sudakov variable defined earlier, kf = kfi + kf2 the transverse 
momentum of the massive parent gluon, z the fraction oi/a and q the relative "transverse 
angle" of the pair. Thus we have: 

ai = za, a2 = {l-z)a, kt = hi + kt2, ^'^■^^^ 

The above can also be straightforwardly expressed in terms of the mass of the parent 
gluon = z(l — z)q'^ so that: 

dT2{k„k2) = n = dm'^^dz^ , (2.12) 

a,- vr an 27r 

i=l * 

where (p is the angle between kt and q. 

We note that the quantity is invariant under boosts in the longitudinal direction 
implying that it is independent of the variable q, which allows us to perform the a integral 
explicitly. Defining the result of the a integral by a trigger function 0, we can write the 
inclusive piece (first two terms on the right-hand side of eq. (2^)) as: 

Cf f d^kt 



{6pty = — {«s(0) + 4nx{k!)} no{k^) + 



with the trigger function being: 

noikf + m^) = [' —5pt (a, /? = (fc? + m^)/{sa)) , (2.14) 

J{k1+rrfi)/s " 



-8- 



and similarly for J7o(A;j) which just corresponds to real massless gluon emission, and the 
superfix i pertains to the inclusive piece. In writing the above we have made use of the 
fact that 5pt depends on the a and /? Sudakov variables and that (3 = {k^ + m?)/{as) in 
the massive case and k^/{as) for massless real emission. 

For the non-inclusive contribution a trigger function can be similarly defined so 
that it reads: 

i,„r = ^c,J{^)\A/^^n„., (2.15) 

with 

fi„i= / ^(Sp^(ki,k2)-dpt{ki + k2)). (2.16) 

J{k1+m^)/s a 

We shall first address the terms involving the trigger functions for massive and 
massless single gluon emission and then focus on the non-inclusive term. 



3. Dispersive calculations 



Here we shall show how the leading non-perturbative terms may be extracted using the 
dispersive representation of the running coupling [l^. We start with the inclusive treat- 
ment, compute the relevant trigger functions and then move to the non- inclusive gluon 
decay term. 



3.1 Naive result and inclusive correction 

Let us consider the terms in eq. ( ^ ) involving ^^ik'^) and nQ{k1 + m'^). As was discussed 
in refs. 21] the integral of the decay matrix element contains a term proportional 



to /So as well as a singular term where the singularity is associated to the collinear gluon 
splitting. The collinear singularity cancels against a similar one in the function x leaving a 
finite remainder. Also the term proportional to /3o combines with the as(0) term to build up 
the running coupling in the CMW scheme p3], which shall always be implied henceforth. 



To be more precise, after integrating over the decay phase-space dT2, eq. (2.13) can be 
written as a sum of two finite terms: 



{5ptY = {sptf + isptr' . 

Let us first describe the term {6pt)^ which reads [21|: 



= / (im^acfr(m^) 



vr 



Cf 



+ m 



dm? Jq k^ + 



"^^^ no {k^ + m') 



dm^ 



m^ 



acsim'^)0,o{m'^). 



(3.1) 



(3.2) 



In writing the above we have, as in refs. [18, 21], made use of the ultraviolet convergence 
of the integrals to extend the upper limits of the integration to infinity, since we are 
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concerned here with the infrared region kt,m Aqcd • We also wrote ~ s as the upper 
hmit of the kt integral, which is redundant here since our final results stem from the low kt 
infrared regime. The above result involves invoking the effective coupling which is related 
to as via a dispersive integral [p^]: 

— dm - — „ „ „ , (3.3) 



so that one has: 



To arrive at the last line of eq. (|3.2|) we substituted the effective coupling instead of ag 
using eq. (3^) and performed an integration by parts. The term corresponding to {6pt)^ 



is the one that is produced by ignoring all two-loop effects except those included via the 
running coupling. It thus corresponds to a naive single-gluon treatment and is thus referred 
to as the naive term. This is the term that will be related to the one-gluon emission study 
in ref. [14|, where a running coupling as{kt) was introduced by hand into the calculation^. 



The term {SptY^'^ is a leftover from the cancellation of singularities between the function 
X and the collinear divergent piece of the result for the integral of |20, 21 1, which 



following refs. [pO| , 21] we refer to as the inclusive correction to the naive result: 



{Sp,)-=^-^[ ^^^^^d [ ^^lnMM±^ff-(fc? + m^), 



/3o J rn? 4tt d In J k"^ + 
O^"'^ = floik^ + m^) - Qoikf). (3.5) 

The inclusive correction merely arises out of the difference between the massive gluon and 
massless gluon trigger functions which we denote by Q'"'^, and is the first sign of two- loop 
effects not included simply via the use of a running coupling. 

3.1.1 Naive result 

In order to focus on the non-perturbative contribution alone we can split = a^"^ + Sag, 
with 6as a non-perturbative modification to the perturbative definition of as which results 
in the corresponding modification to the effective coupling Jaeff- 

We can also take the small-m^ limit of the trigger function i^Q^m?) so that the non- 
perturbative (NP) contribution reads: 



{^Pt)t^p = — [ ^^^^6ac{i{rn^)^o{m'^). (3.6) 



The non-perturbative modification 6a must vanish in the perturbative region so that the 
coupling Q!cfT matches on to the perturbative coupling. Hence the integral above, represent- 
ing the non-perturbative contribution, has support only over a limited range of up to 



^Strictly speaking the running coupling used in ref. had as argument the transverse momentum with 
respect to the dipole axis in the rest frame of each emitting dipole. For generating the 1/R collinear term 
it suffices to replace this with kt the transverse momentum with respect to the emitting hard leg. 
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some arbitrarily chosen perturbative matching scale ///. The /U/-dependence of the result 
is then cancelled up to 0{a'^) when combining the non-perturbative result with a pertur- 
bative estimate at the same order. This issue, mentioned for the sake of completeness, is 
irrelevant to the rest of our discussion and shall henceforth not be mentioned again. We 
shall also avoid writing the subscript NP since we always refer to the non-perturbative 
contributions from now on. 

Let us compute the naive trigger function ^o{k^ + m^) in order to evaluate eq. (^). 
This is given by: 

no{k^ + m^)= [ —6pt{k), (3.7) 

J{k1+m?)/s « 

where = fci + A;2 is the parent (massive) gluon four-momentum. Using the results of the 



previous section eqs. (|2.3|) to (|2.5| ) we have: 



5pt{k) = Spj + 5p^ 

= -^/3ein(fc)-^aGout(fc), (3.8) 

where Bout and ©in are step functions indicating whether the gluon is outside or inside the 
jet (i.e. if it is recombined with the hard parton or not). 

In all the jet algorithms that we consider, the clustering of a single particle (e.g. a 
soft gluon) to the triggered hard parton is based on the distance criterion d'rf + 54)^ , with 
5r] and dcj) the separation in rapidity and azimuth between the particle considered and the 
hard parton. At small angles (i.e in the collinear limit relevant to our results) this is just 
9"^, where 6 is the small angle between the soft gluon and the triggered hard parton. Thus 
at the single gluon level the soft gluon is combined with the hard parton if 9"^ < i?^, and is 
not combined if the converse is true. Moreover using the easily-derived relation^ 9"^ = 4ij3/a 
we have: 



eout(A:) = e (4/3/a -R^) = Q [2^^^ 



a , (3.9) 



R^s 

where we used a/3 = (k\ + rn?)/s, and conversely for <d\^{k). 

Let us first consider the ©out = 1 case or the gluon not recombined with the jet. The 
contribution to the trigger function is: 

Q„« + ™^) = -f/' ^„ef2V;|pf-„V (3.10) 



Performing the above integral we arrive at: 



(3.11) 



^Strictly speaking this relation is true for massless gluon clustering. In the nai've calculation we are 
however free to ignore the gluon mass in the definition of the clustering angle, since precisely the same 
quantity is subtracted in the inclusive and non-inclusive terms at the end. This freedom to define the role 
of mass in the naive calculation was also exploited in refs. ji^, ^ to obtain a trigger function oc ^Jkl + m?. 
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where we retained only the leading 0{kt/R) term we wish to compute, neglecting terms 
of order fcj that arise from the lower limit of the integral. Thus the 1/R behaviour is 
associated to quasi-collinear soft emission such that 1 » a » /3, where 1 » a corresponds 
to soft emission, while a ^ /? is the collinear regime. We note that the convergence of the 
a integral can be used to extend the integral to infinity, while the irrelevant lower limit can 
be set to zero as in refs. [20, |2T|, J da/a ^ da/ a. Lastly we just mention that the term 



with the soft gluon recombined with the jet, 5pJ , does not contribute to the 1/R result. 
The physical reason for this is straightforward: the contribution to the jet 6pt from this 
term is proportional to [3 or the jet mass, which vanishes in the collinear limit responsible 
for the 1/R behaviour, i.e. is collinear-safe. Thus only gluons not recombined to the hard 
jet give us a hadronisation correction varying as 1/R. 

Now we substitute eq. ( |3.11| ) into the result eq. ( |3.2| ) to obtain the naive result: 

vn 1 Cp f°° dm? ^ ,9, , , 

{Sptf = —B— -^6a,s{m^)m. (3.12) 
n IT Jq m 

This result can be written as 

24i 

(5pt)° = — ^, (3.13) 
where we introduced the non-perturbative coupling moment Ai: 

Cf f°° dm^ . , 9^ . ^ 

^1 = ^/ -^m6a,sim^). (3.14) 





The coefficient 2Ai of the —1/R term can be compared directly to the naive result 
for e+e~ event-shape variables pO| , 21 1. It is precisely one-half of the result for the thrust 
mean-value exactly as was reported in IQ also within a naive single gluon approach. 

3.1.2 Inclusive correction 

Next we deal with the inclusive correction term. Here the trigger function is just the 
difference between the massive and massless gluon phase-space: 



= f 2^^ - « U ^ f f ^ - a) 

Jo a [ R^s i 2 Jo a \R^s J 



2 



^[\lk^t+rr?-h]. (3.15) 



The inclusive trigger function is precisely of the form one meets in refs. |2C, 21 1, with 
the coefficient of —1/R, being identical to the corresponding quantity for event shape- 
variables, enabling us to simply use the results derived for event shapes. Inserting the 
trigger function into eq. ( |3.5D and introducing Sacsirr?) as before, one gets: 

{6pty'''={6pt)\^n, (3.16) 



where Tin = 3.299C^//3o is the value quoted in refs. 21|. The factor ri„ is the same 



as for the thrust and other event shapes so adding the inclusive correction to the naive 
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result, the result for the jet 5pt is still one-half of the corresponding result for the thrust, 
i.e. universality is maintained at the inclusive level. 

Having computed the naive and inclusive pieces we are left with the non-inclusive 
correction. It is here that the jet 6pt result will no longer receive the same correction as 
event shapes and universality of the Milan factor will be broken (except for the case of the 
anti-/cj jet algorithm). 

3.2 Non-inclusive correction 

For the non-inclusive term one can follow a similar treatment to the inclusive case but with 
the appropriate non-inclusive trigger function. Then eq. ( |2.9D becomes, as for the case of 
event-shape variables |20, 



/A \™ /•°°rfm^ aeff(m2) d f f^^ dc/, f\ f . ^ 



Po Jq m? Att dlnrn^ \Jq 2tt Jq J 2! 

We shall focus our attention on the term in parenthesis that yields the linear behaviour in 
m corresponding to the leading non-perturbative (hadronisation) correction. In order to 
compute this we need to evaluate the non-inclusive trigger function: 

r°° dn 

/ —{6pt{k^,k2)-6pt{k, + k2)). (3.18) 
Jo a 

The term 5pt{ki, /C2) represents the contribution from two parton emission, which depends 
on the details of whether the offspring partons end up inside or outside the jet after the 
application of the jet algorithm. For most jet algorithms currently in use, this involves 
more than just working out whether the gluons are within a distance R of the hard parton 
(an exception is the anti-/cj algorithm), so we shall need to consider the action of the 
algorithms in some detail. Due to the length of this discussion we shall devote the next 
section entirely to it. 



4. The non-inclusive trigger function 

We shall start by considering the kt algorithm. Let us briefly remind the reader of how 
this works. 

First one computes the distances d^j between all pairs of objects i and j, dij = 
min ^K^j^,K^j^ Uij, with yij = 5rf + 5<j)^ , where Sr] and 5<l) refer to the separation of i 
and j in rapidity and azimuth measured with respect to the beam and we used Kt to de- 
note the transverse momentum with respect to the beam direction (while reserving kt for 
the transverse momentum with respect to the hard emitting parton). In the small angle 
approximation appropriate to our work we just have 5rf + 5(p'^ 9fp where 6ij is the angle 
between i and j. One also computes the distances dis between each object and the beam 
defined as ref j-R^, with R the selected radius parameter. 

Then if amongst the various distances the smallest is a diB, object z is a jet and is 
removed from the list of objects to be clustered. If the smallest distance is a d^ then 
i and j are combined (merged) into a single object and the procedure is then iterated 
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until all objects have been removed. The recombination scheme we use here is addition 
of four-momenta so that the object resulting from combining i and j has four-momentum 
Pi+Pj- 

We note from eqs. (|2.3| ) to 
of ki and k2, 

Spt{ki, k2) = 5pt{ki) + 5pt{k2) 

= 6pt{ki)Ein{ki) + 6pt {ki)Eont{ki) + {ki ^ /ca) , (4.1) 

where we expressed 6pt in terms of the contributions from ki and k2 ■ We also distinguished 
the cases when the offspring partons are recombined into and fly outside the jets as the 
contributions are different in either case: 

5pt{k^) = -^P^, 

6pt{h) = -^ai, (4.2) 

where pf expresses the contribution from a recombined parton and for a non-recombined 
emission. The conditions Hin,out denote the constraints for partons to be in and outside 
the jets and they are no longer simple step functions as for the single gluon case, and 
must be obtained by applying the algorithm in full including the possibility of soft partons 
clustering to each other. 

4.1 6pt in the kt algorithm 

Let us first apply the kt algorithm to work out the 6pt{ki,k2) in various situations. We 
have to consider the soft partons ki , k2 and their distances from the triggered hard parton 
"jet" which we denote by j. To this end let us define the step functions = © [R^ — O^i^) , 
so that if two partons have an angular separation less than R the step function is unity 
and else it is zero^. Then we can divide the full phase-space up into definite regions (as 
indicated by table |l]) which account for all possible cases for the angular separations. For 
each row of the table (i.e each of the possible combinations of values of the 0^6, it is possible 
to obtain a value for the 5pt{ki, k2). 

In table ID we have written 5pf~^{ki^ k2) = 5pf {ki)+5pf {k2) to indicate the contribution 
when both ki and k2 are eventually combined into the jet, while dp^ {ki, fca) = 8pf{ki) + 
6p^{k2) indicates a situation where ki is recombined and k2 is not. Let us show how 
the entries in this table are obtained by considering an example and discussing the more 
interesting scenarios. 

The scenario: ©12 = 1, ©ji = 1, Qj2 = 0. We consider the case described in the second 

^Recall that we are looking at the small-/? behaviour R <^ 1, which entitles us to replace St]^ + 54>^ by 
the small angle approximation 9^. Neglected terms produce only finite R corrections. 



2.5) that the 6pt{ki,k2) is additive over the contributions 
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0,1 


0,2 


Spt{ki, k2) 


1 


1 


1 


Sp++(ki,k2) 


1 


1 





Discussed below 


1 





1 


Discussed below 


1 








Discussed below 





1 


1 


Spt^{ki,k2) 





1 





^pt~(.h,k2) 








1 
















Table 1: Values of Spt depending on the angular separations of partons. The three entries 
"Discussed below" need special treatment and thus are left for discussion in the text. 

row of the above table. The distances we need to discuss are the following: 



diB 


- k2 K'2 

— Kf irL , 


d2B 


— fv(2-"' ' 


di2 


■ / 2 2 \/)2 


dij 


— '^t,l"lj ) 


d2j 


— '^t,2t'2j • 



In this situation (given the values of the 6 functions) it is not possible for any of diB, 
d2B, or d2j to be the smallest. Thus either di2 or dij is the smallest. 

If di2 is the smallest (i.e. smaller than dij) then ki and k2 are clustered to each other 
first. Two sub-scenarios arise: the resultant soft jet may fall either within or outside 
the grasp of the hard parton. If the soft jet is recombined with the hard parton the 
resulting contribution is Spf'^, else it is Spl . 

Alternatively if dij is the smallest distance particle ki is combined with the jet first. 
Then particle k2 forms a separate jet, i.e. is removed as it is closer to the beam than 
it is to the hard jet. The contribution is then Sp^ {ki,k2)- We can summarise the 
overall scenario just discussed by the result: 

Case di2 < dijt 

Spt{ki,k2) = e{dij - du) [e{R^ - olj)5pt+ + 6(-i?2 + elj)Sp;-] , (4.3) 

where 6kj is the angle of the soft jet (parent gluon) k = ki + k2 with the hard 
parton. 

Case di2 > dij: In this case particle 1 gets clustered to the jet and 2 is left outside. 
Thus: 

Sptih, k2) = e{di2 - dx,)6pt-{kiM)- (4.4) 
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The scenario: Q12 = 1? ©ji = 0, Qj2 = 1. This scenario (third row of the table above) 
is the same as the one just discussed but with ki and k2 interchanged. 

The scenario: B12 = 1, ©ji = , Gj2 = 0. Here the smallest quantity is di2 and the soft 
partons are clustered to each other first. The resulting soft jet may or may not be 
recombined with the hard parton. Thus we have in this region: 

6pt{kuk2) = e{R^ - e%)6pt^ + e(-i?2 + e%)5pr , (4.5) 

where the step functions specify whether the soft parton jet is within an angular 
distance R of the hard parton or not. 

All the other scenarios follow straightforwardly and the values of the 5pt are as indi- 
cated in the table. For instance if one considers the bottom row of the table where all the 
Q functions vanish (partons always separated by more than R) the diB are always smallest, 
and the partons ki and k2 are never combined in to the triggered jet, so the result is 6p^ . 

Now we are in position to write down the expression for 6pt at two-parton level. We 
express it as follows: 



Spt{ki,k2) = Spf^{ki 

+ Spt+{k^ 
+ Spt+{k^ 

+ spi'iki 



k2)ein{kl)@in{k2) + 

k2)Qont{ki)eontik2) [1 " ei2(A;i, A;2) + ei2(fci, A;2)eout(A;)] + 
A;2)©out(A:i)eout(A:2)ei2(A:i, A:2)ein(A:) + 

k2)einikl)eoutik2)@12{ki,k2)e{dij - di2)ein{k) + (fci ^ k2) + 
k2)Qm{kl)@out{k2)&12{kl, k2)e{dij - di2)@out{k) + {ki ^ k2) + 

k2)Qin{ki)<doAk2)Qi2{ki, k2)e{di2 - dij) + {ki ^ /C2) + 
A;2)©in(A:i)eout(A:2)(l - ©12) + (^i ^ ^2). (4.6) 



In the above ©in(^j) = &{R^ — Ofj) and ©out = 1 — ©in pertains to whether the parton ki 



is within or outside an angular separation R of the hard parton, while ©12 = @{R^ — 0^2) 
is the constraint for the soft partons to have an angular separation less than R. Eq. ( [4.6| ) 
merely summarises the contents of the table, for instance the sum of entries in the first 
and fifth rows is just the first line of the above equation. For consistency one can check 
by removing all the 5pt terms in the above equation and adding all the theta functions as 
they appear one should get 1. The check is positive (meaning all possible scenarios have 
been discussed without double counting). 

We can also write eq. ( [4.6| ) in terms of the individual contributions of ki and /c2 since 
the 5pt is additive over the contributions of individual partons: 

Spt{ki,k2) = 5pt{ki) + 5pt{k2) = 5p~[ {ki)'~ont{ki) + dpf {ki)'~;r,{ki) + {ki ^ ^2), (4.7) 

where we separated the contributions from each offspring parton into its + and — compo- 
nents according to whether the parton is recombined with the jet or not. 

Next we note that for the 1/R term one is only interested in partons not recombined to 
the jet. Partons recombined with the jet contribute to the 6pt via 6pf{ki) oc Integrating 
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such contributions over the da /a phase-space produces only terms regular in R, as for the 
single gluon case. Thus for computing the 1/R term the 6pt{ki, ^2) can just be expressed 
as: 

5pt{ki,k2) = 5p^ {ki)Eout{ki) + Sp^ {k2)Eout{k2) ■ (4.8) 
Thus from eq. ( [4.6| ) focussing on all terms where ki is not in the jet we can identify: 

Sout(A;i) = eout(A:i)eout(A;2) [1 - Qi2{kuk2) + ei2(A:i, fc2)eout(fc)] + 
+ ein(A:i)eout(A;2)ei2(A:i, A;2)G(dij - (ii2)eout(A:) + 
+ ein(A:2)eout(A;i)ei2(A:i,/c2)G(d2i - c^i2)eout(A:) + 
+ ein(A;2)eout(A;i)ei2(/ci, A:2)0(di2 - d2j) + 

+ ein(A:2)eout(A:i)(l-ei2), (4.9) 
which can be simplified after some algebra to: 

Hout(A:i) = eout(A:i) [1 - eout(A:2)ei2(A:i, A;2)ein(A:)] + 

+ ein(/ci)eout(fc2)ei2(/ci, /C2)e(dij - (ii2)eout(A:) + 

- eout(A:i)ein(A:2)ei2(A:i, A;2)e(d2i - di2)@in{k), (4.10) 

which represents all possible ways for ki to be outside the jet after application of the kt 
algorithm. A similar expression holds for k2. 

Repeating the above procedure for the anti-Zct algorithm yields a much simpler result 
due to the absence of soft parton self-clustering. Thus in the anti-kt case a given parton ki 
is always outside the jet if it is separated by more than R in angle from the hard parton: 

Hout(fe^) = Gout(A:^). (4.11) 

4.2 The a integral 

Having worked out the relevant piece of the 5pt{ki, ^2) for the kt algorithm we now proceed 
to carry out the integral over a involved in the definition of the trigger function eq. (|2.16| ) . 
This is given by: 

/dcy, 
— {6pt {ki)Eont{kl) + Spt {k2)Eont{k2) - Spt {k)Qont{k)) 

/s f d(x 

= — TT i — (aiHout(fci) + a2^ont{k2) - aQoutik)) 
2 J a 

= — — / da (2;Hout(fci) + (1 - z)Eontik2) - Gout(fc)) 

= ^1 + ^2- ^3, (4.12) 



where 



^1 = -\z da[Qout{ki) [1 - eout(A:2)ei2(fci, A;2)Gi„(A:)] + 



+ @in{ki)eont{k2)ei2ikl,k2)&{dij - di2)@out{k) + 

- Gout (fei ) Gin (fe) 912 (A:i,A;2) 9(^2, - di2)@in{k)], (4.13) 
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with eg. 0out(^i) = 0(^ij — R^)- Note that in writing the above we have used 5p~[{ki) = 
—aiy/s/2, with ai = za and 02 = (1 — z)a. We shall give expressions for and $^3 in 
due course. 

In order to proceed we shall also need to express the angles between partons in terms 
of the Sudakov variables. The angle between the soft parton ki and the hard parton j is 
just (assuming ai ^ Pi, which is the collinear approximation) Ofj = 4/3i/ai = 4A;^]^/(saf ). 
Moreover since ai = za, Ofj = Aqf/{sa'^), where qi = kti/z. Similarly we have = 
4g|/(sa^) with q2 = kt2/{l — z). The angle between the soft jet formed by clustering ki 
and k2 and the hard parton is given (again in the collinear limit) by 9\- = 4:kf/{sa'^), and 
finally the angle between ki and k2 is just Of 2 = 4(7^/ (sa^), where q^ = {kti / z — kt2 / — z))^ ■ 
In writing these angles we have used the small angle approximation throughout, specifically 
(1 — cos 9) ~ O"^ /2, and also considered ai ^ /?, as we are always concerned with partons 
in the region collinear to the emitting hard parton. Terms neglected in this approximation 
contribute to corrections regular in R. 

Next we re-scale ^Rajl — > a so as to straightforwardly extract the XjR dependence 
that concerns us. Hence we write: 

= -\-z [ da{e {qi - a) [1 - e {q2 - a) Q {a - q) Q {a - h)] + 
^ Jo 

+6 (a - qi) @ {q2 - a)@{a- q)Q {kt - a) Q {Kt,iQi ~ [i^t,!^ '^4,2] 9^) + 
-9 {qi - q) e (a - 92) e (q - g) G (a - h) & {t4^2l2 - min ^^,2] Q^) } 

= ~ (min[gi,g2] - max [g, fcj]) 9 (min [gi, ^2] -max[g,A;t]) + 

+ (min [^2, h] - max [qi,q\) 9 (min [q2, h] - max [qi,q]) 9 [qf - q^miu [l, k^^2/'4,i\) + 

- {qi -ui&y.[q2,q,kt])Q{qi - max [^2, 9, A;*]) 9 [ql - q^mm. [k^^i/k? 2, l] ) }, (4.14) 

where we used the fact that the distances dij between partons involve the transverse mo- 
menta K with respect to the beam direction which must be distinguished from kt, the 
transverse momentum with respect to the jet direction, and we performed the trivial a 
integral extracting the crucial 1/R dependence we seek. 

Note that the transverse momenta with respect to the beam are essentially the lon- 
gitudinal momenta with respect to the jet (recall that we are looking at jet production 
at ninety degrees to the beam as well as parton emission very close to the triggered hard 
parton, Ui ^ (3i), so we can to our accuracy, just use Kt^i ~ y^/2ai = y^/2za and 
likewise for Kt,2 with z ^ 1 — z. 

Thus we can express eq. ( [4.14 ) as: 

Ql = -^^{91 - (min [ft: ^2] - max [q, kt]) + 

-1-^ (min [q2, h] - max [qi,q]) x 9 (gi - gmin [1, (1 - z)/z\) + 

-■^ (gi - max[g2,g, fct]) x 9 (92 - gmin[2;/(l - z),l]) |, (4.15) 

where we defined the function ^{x) = xQ{x). Here ^'(2;) = x if x > and ^'(x) = 
otherwise. 
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A similar result is obtained for the dpt{k2) contribution: 

O2 = ~ - ^ (min tell 92] - max [q, h]) + 

(min [qi,kt] - max [q2,q]) x & {q2 - qmm [l,z/{l - z)]) + 

- max [qi,q, h]) x e {qi - qmin [(1 - z)/z, 1]) }. (4.16) 

The subtracted f^s contribution is simply the subtraction of the naive trigger function 
computed before: 



1 

R 



(4.17) 



The overall result for the non-inclusive trigger function is then given by combining the 
ill, ^^2 and terms and takes the form: 



j^(^hi+kt2- ^k^+m^ + f{qi,q2,z)^ , (4.18) 



which we wrote in terms of the ku for sake of easier comparison with the corresponding 
expression eq. (3.31) in ref. pl[ . The first three terms of the result correspond to the 
non- inclusive trigger function computed for event shapes while the extra term f{qi,q2,z) 
will break the universality of the Milan factor. 

To express the result in terms familiar from refs. pO| , ^] we introduce ui and U2, where 
= Qi/Q- Defining Q,{ui,U2) so that 0,ni{qi,q2) = —q/R^{ui,U2) we arrive at: 



J7(iii, M2) = zui + (1 - z)u2 - J zu{ + (1 - z)ul + 



^' ( min[Mi, U2] — max 



1, \J Zul + (1 — Z)U2 — z{l — z) 



+ 



+ z-^ ^' { min 



U2, \J zul + (1 ~ ^)'^2 ~ ~ ^) 

x9 (ui — min [1, {I — z)/z\) + 



max [til, 1] I X 



^' I ui — max 



U2, 1, \J ZU\ + (1 - Z)u1 - Z{1 - z) 



xQ{u2- min [z/{l - z), 1]) \ + 



+ (1 - z)< ^ ( min 



Ui, Zu\ + (1 — Z)U2 — z{l — z) 

xQ {u2 — min [z/{l - z), 1]) + 
— ^ I U2 — max 



max [u2, 1] I X 



til, 1, ^J zul + (1 ~ ^)'^2 ~ -^i^ ~ A 



x6 (ni — min [1, (1 — z)/^]) 



(4.19) 



where we used (fcj + m?)/q^ = zu\ + (1 — z)u\, which follows from the fact that m? 
z{l - z)q^. 
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Once again the first line of the above result corresponds to the Milan factor for event- 
shape variables, and is the same trigger function as encountered in refs. ||2^, 21 1, while the 
remainder will constitute the correction to it from kt clustering. 

The result corresponding to the anti-A;^ algorithm is also just the first line of the above 
result ( ^.191 ) since there we have: 

^i = -^z [ — aeiqi-a), (4.20) 
R J a 

which gives —q/R (zui). Similarly the 0,2 term is the same with z ^ 1 — z and the naive 
piece is of course common to both algorithms. 

Thus in the anti-A;t algorithm after combining all the pieces we obtain that the result 



will be identical to that for event-shape variables, computed in refs. [20, 21]. 

Having derived the trigger function in the kt algorithm we now proceed to perform the 
integration over the remaining variables, including also the decay matrix element as 



required by eq. ( p. 17 ). The details of the integration procedure will be consigned to the 



appendix and in the next section we shall quote the result and combine it with the naive 
and inclusive terms. 

5. Results and conclusions 

In this section we will mention the result obtained for the non-inclusive correction using 
the trigger function of the previous section. We shall then combine the result with the 
naive and inclusive pieces computed before. 

After computing the trigger function described in the previous section one inserts 
the result into eq. ( 3.17| ) to obtain the result for {6pt)^\ In order to do so we need to 



integrate over the parton decay phase-space including the squared matrix element 
describing gluon decay. The details of this integration are mentioned in the appendix and 



here we just quote the result analogous to eq. (4.8) of |21]: 

{dptr = {6pt)'rm, (5.1) 

where for the kt algorithm we have a result for r„j different from that for the event-shape 
variables studied in refs. |^, 

rni = J-(-2.145Ca + OMOCa - 0.097n/) 

= -^(-3.071Ca - 0.194n/), (5.2) 

Po 

where the results in the first line correspond to the "soft gluon", "hard gluon" and "hard 



quark" contributions in the terminology of refs. [20, 21]. 

The corresponding result for the anti-/cj algorithm coincides with that obtained for 
event shapes [^, 

r„i = A(-1.227Ca + 0.365Ca - 0.052n/). (5.3) 

PO 
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Combining the result obtained for r„j with that for the inclusive piece (eq. (l3l^ )) and 
the naive result (eq. (|3l3|) ) one obtains the full result for jets defined in the kt algorithm: 

-2Ai 

{6pt) = ^ (1 + Tin + rni) 

= ^(1-01), (5.4) 
for Uf = 3 (for Uf = the rescaling factor is 1.06). The factor multiplying the naive 



result is thus found to be 1.01 compared to 1.49 for event shapes pQ , 21| and the anti-Zc^ 
algorithm. Here we took n/ = 3 since we are dealing with the non-perturbative region but 
in the corresponding perturbative estimates one should of course take Uf = 5. To obtain 
the results for a gluon jet one can just replace Cp in Ai, by Ca- 

Thus compared to the thrust 1/Q correction which has a Milan factor of 1.49 the 
1/R term of the jet hadronisation for quark jets in the kt algorithm has a coefficient that 
is 0.5 X 1.01/1.49 ~ 0.34 times the coefficient of 1/Q for e~^e~ thrust instead of 0.5 as 



observed at the single gluon level in [14|. For the anti-kt algorithm the result is still one- 
half of that for the thrust after inclusion of the Milan factor 1.49 in both cases. The 
somewhat smaller hadronisation correction predicted here for the kt algorithm compared 



to the anti-A;^ algorithm was also observed in the Monte-Carlo studies of Ref. |14|. 



To summarise, we began this paper by noting that the single-gluon calculations per- 



formed in ref. |14| are inadequate in terms of determining the size of the hadronisation 
correction to jet pt relative to the known 1/Q corrections for event shapes. The reason 
for this inadequacy was merely the fact that a two-loop calculation had been shown to 
be necessary in refs. ||2^, ^ to unambiguously determine the size of 1/Q hadronisation 
corrections to event shapes. 

Our aim in this paper was to carry out a similar calculation for the jet pt so as to finally 
put it on the same footing as event shape variables, which makes the comparison feasible. 
We carried out such a calculation and pointed out that as for event-shape variables one 
can write the result as a two-loop rescaling factor multiplying the naive single-gluon result. 
Unlike event-shape variables however, for jets defined in the kt algorithm the rescaling 
factor is not the universal Milan factor 1.49 turning out instead to be 1.01. We confirmed 
the expectation of ref. |^] that for jets defined in the anti-fc^ algorithm the rescaling factor 
is 1.49, as for event-shape variables. 

We also note that we have presented here the calculation of correlated two-parton 
emission and neglected higher correlations involving three or more partons which would 
come in beyond our two-loop accuracy. However the two-loop calculation we performed 
here corresponds to the accuracy needed to fix the scale of the coupling or equivalently 
to fix a value for Aqcd which directly controls the size of the non-perturbative contri- 
butions. Hence while one expects the eff^ects of the two-loop calculation to be crucial, 
higher-loop calculations can be expected to contribute at most corrections of relative order 
Os/tt which could change the value of the two-loop factors at most at the twenty percent 



level [20|. Interestingly Sterman and Lee [25| have recently shown within the context of 



soft-collinear effective theory that for observables such as event-shape variables that are 
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linear in emitted parton transverse momenta, one can generally demonstrate universality 
of power corrections without resorting to fixed-order perturbative calculations beyond the 
naive massive gluon results. It would be of interest to revisit their arguments in the context 
of the observable studied in the current paper. 

To conclude we point out that the potential experimental studies mentioned in |14| 
should of course take into account the result of the calculation performed here for jets 
defined in the kt and anti-fct algorithms. The calculation for jets defined in the SISCone 



and the Cambridge/ Aachen algorithms is work in progress |2£] 
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A. Numerical result for the non-inclusive correction 

In this section we evaluate the integral relevant to computing the non-inclusive correction. 
Following ref. [^], one can write eq. ( |3.17| ) as: 



with: 



7r/3o 



dz 



du- 



,2 J I 



duA 



--U1U2 



1 



zu{ + (1 — z)u\ 



(A.l) 



(A.2) 



Here we have defined u± = U2 and used eq. (4.14) of ref. [^]. The squared matrix 
element is given in eqs. (2.8) and (2.9) of the same paper pO|]. 

By symmetry we just consider the region of the phase-space of integration defined by 
u_ > (i.e. ui > U2) and multiply the result at the end by a factor 2. This way the trigger 
function (eq. ( [4.19| )) takes the simpler form {ui > ^2): 

(l{ui,U2) = 



ZUl + (1 - Z)U2 - \J Zu\ + {I - z) 

—'^ I U2 — max 



ul + 



z^ I ui — max 



1, Y Zu\ + (1 — z)u'2 — z{l — z) 
U2, 1, \J ZU\ + (1 — Z)U2 — Z{\ — z) 



+ 



6(^2 -min [z/{l - z),l\) + 



+{l - z)^ Jzul + {1- z)u\ - z(\ - z) - max [u2, 1] 6 (^2 - min [z/{l - z), 1]) , 



(A.3) 



where we made use of the easily-derived inequality U2 < \/ zu'f -|- (1 — z)u2 < ui, for 
< z < 1, and hence y^i?I^~-i-~(r^^^)ti^"^^^(r^^^ < ui. 

To further simplify the trigger function one needs to consider in detail the relative 
size of the various components of the functions ^ above. To this end we proceed in the 
following way. 
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A.l Simplification of the trigger function 



We consider all possible permutations of the relative size of ^zu{ + (1 — Z)U2 — z{l — z), 
ui, U2, and 1. Thus we discuss the following scenarios: 



1 < U2 < Zul + (1 — Z)U2 — z{l — z) < Ml or 1 < Zul + (1 — z)u2 — z{l — z) < 

U2 < ui. We call the corresponding region of the phase-space "A", thus: 



^a{ui,U2) = y + (1 — z)u2 — z{l — z) — \ zu\ + (1 — z)u 



(A.4) 



■\/ zu'f + (1 — z)u'2 — z{l — z) < 1 < U2 < ui. We label the corresponding region 
Here we have: 



U2) = 1 - a/ Zu{ + (1 - Z)U2 ■ 



(A.5) 



U2 < 1 < a/ zu\ + (1 — z)u2 ~ ^(1 ~ z) < ''^i- We name this region "B". In this case 
we have: 



(Ib{ui,U2) = ZUi + (1 - Z)U2 - J zuf + (1 - 2)^2 + 



+ 



Zuf + (1 — Z)U2 — Z{1 — z) — ZUl — {1 — z) 



Q{u2- z/{l- z))e{l-2z). 

(A.6) 



U2 < \J Zul + (1 — Z)U2 — Z{1 — z) < 1 < Ml or Y^Zuf + (1 — Z)U2 - z{l — z) < U2 < 

1 < ui. We call this region "C", hence: 



Q.c{ui,U2) = ZUl + (1 - z)u2 - y zuf + (1 - z)ul+ 

- z (ui - 1) e {U2 - z/{l - z)) 6(1 - 2z). (A.7) 



U2 < \/ Zu\ + (1 — z)u\ — ^(1 — z) < Ui < 1 or \J Zu\ + (1 - z)u\ — z{\ — z) < U2 < 

ui < 1. We name this region "D", where: 



QDiui,U2) = ZUl + (1 - Z)U2 - J zul + (1 ~ ^)'^2 ■ 



(Ai 



We summarise the above results for the trigger function in table Q for which we define 
Clm = ZUl + (1 — z)u2 — \/ zu'f + (1 — z)u2, which is just the non-inclusive trigger function 
for event shapes and the anti-A;^ algorithm. We also show in fig. || the {ui,U2) phase-space 
and the relevant regions of integration. We note that the overall region of integration is 



u+ > 1 and —1 < n_ < 1 [20|, but we are only considering the region < u_ < 1 as we 
stated before. 
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Region 


Q{ui,U2) 


A 


J zu\ + (1 — z\uk — z(l — z\ — J zu\ + (1 — z^uk 

= + \/ -2^1 + (1 — Z)'^ ~ -2(1 — -2) — ZU\ — (1 — z)u2 


A 


1 — zuf + (1 ~ z)Xi^ = Vlra + 1 — ZU\ — (1 — z)u2 


B 


+ V Zuf + (1 — Z)U2 — z{l — z) — ZUl — — z) X 

xe[u2 - z/{i - z)]e(i/2 - Z) 


C 


n^-z {ui - 1) @[u2 - z/{i - z)]e{i/2 - z) 


D 





Table 2: The trigger function in different regions of the phase-space. See fig. 




Figure 1: The (ui,U2) phase-space of integration for eq. (A2). See table || for the trigger 



functions f2 in various regions. Here the elhpse (curved hnc) defines the boundary between regions 
governed by the relative size of zu\ + (1 — z)u\ — z{l — z) and 1. This ellipse ranges from being a 
straight line parallel to ui axis (when z = 0) or U2 axis (when z = 1) to a circle of radius -\/5/2 
(when z = 1/2). 

A. 2 Convergence issue 

Having simplified the trigger function we now check the convergence of the integral ( |A.2| ). 

First of all we discuss the limit if+ — > 00. The trigger function in this limit (corre- 
sponding to region A) is given by: 



Q.A = \J zul + (1 — z)u\ — z{l — z) — zu\ + (1 — z)u'2 

2z{l - z) 



M-f — ►OO 



u+ + {2z - l)u_)2 + 4z(l - z)(n^ - 1) + {u+ + {2z - l)n_)2 + 4z{l - z)u^_ 

z(l — z) , . s 
-, A.9 
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which clearly vanishes as u+ ^ oo, thus ensuring convergence for reasons cited in |20|. 

Next we discuss the limit z{l — z) ^ and verify the convergence in various regions 
of the (ni,M2) phase-space. Dealing first with region A we notice that when z{l — — > 
the trigger function has the behaviour: 

^^.,1-^-0 41 -z) 

n+ + [2z — 1)U- 

which removes the singularity in the matrix element of the form 1/^(1 — z) (see ref. [pC|]). 
We do not here encounter the problem of the edge of the phase-space as in {u+ — > 1 
and u- ±1) since we are away from these corners of the phase-space in region A, i.e. 
the denominator in eq. ( A.10| ) never approaches zero. 



In region A the trigger function has the following form: 



n^=l-^^zui + {l-z)ul. (A.ll) 

Here we have the condition zuf -|- (1 — z)u2 — z{l — z) < 1, which can be written as 
— fijj < z{l — z)/2 + 0{z^{l — z)"^) when z{l — z) — > 0. Thus the trigger function here 
removes the singularity in the matrix element of the form 1/^(1 — z). 

Next we consider region B. The convergence of the integral involving has been 



discussed in [20|, so we only check the convergence of the integral involving the "correction" 



term 0^'''' = \/ zu\ + (1 — z)u'2 — z{l — z) — zui — {1 — z), which is present if z < 1/2 (i.e. 
we only worry about z ^ singularity). In this region we have the condition zuf -|- (1 — 
z)u2 — z{l — z) > 1, which can be written as < {ui — l)z, where < ui — 1 < 1. 

Thus we deduce that the trigger function removes the singularity l/z in the matrix element 
when z — > 0. 

Similarly for region C the correction term (which is present if z < 1/2) is: —z{ui — 1), 
which clearly removes the singularity 1/z in the matrix element since < -ui — 1 < 1. 

Lastly the integral in region D is obviously convergent since the trigger function there 
is equal to Hence we deduce that the trigger function always removes the l/z{l — z) 
singularity in the matrix element. 

We note that here we do not specifically have the problem of the edge of the phase- 
space which appears in the Milan factor computation ||2^ (except for i}^, which has been 
discussed in the same paper). 

A. 3 Numerical result 



Inserting the results for the trigger function from table ^ into eq. ( [A.2 ), taking into account 



the appropriate phase-space, we are able to numerically compute the result for the non- 
inclusive factor r„j. We present here the numerical result for each region of the phase-space. 
We first compute the result for integration over 0^ which appears in all the regions 



(see table HI). This has already been done [20| and the result is: 

1 



^ , 1.227Ca + 0.365Ca - 0.052n/), (A.12) 
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where after multiplying by the factor 2 (accounting for the fact that we have considered 
just U- > 0) we arrive at the result in [^] . Next we compute the "corrections" region-by- 
region. 

Region D has no corrections, so we begin with region C: the correction term is 
—z {ui — 1), and the phase-space is Q[u2 — z/{l — z)]Q{l/2 — z)@{l — zu\ — (1 — -|- 
z{\ — z))Q{v?L — (m+ — 2)^), with < U- < 1 and 1 < n-|- < 3. The numerical result in this 
case is: 

r^. = J-(_o.262Ca + 0.072CA - O.OOSn/). (A.13) 

The correction term for region B is \/ zu\ -|- (1 — z)u2 ~ ^0- ~ z) — zui — {l — z), and the 
phase-space is e[u2-z/{l-z)]e{l/2-z)e{-l+zuj + {l-z)ul-z{l-z))e{u'i-{u+-2f), 
with < u_ < 1 and 1 < U-|- < 3. The numerical result is: 

r^. = J-(_o.043Ca + O.OUCa - 0.002n/). (A.14) 

For region A the correction is 1 — zui — (1 — z)u2 and the phase-space is 0((u_|_ — 2)^ — 
ui)Q{l- zuf - {l-z)ul + z{l-z)), with 2 < u+ < l + \/2 and < u_ < 1. The numerical 
result is: 

r^. = i-(_o.003CA + O.OOICa - O.OOOn.). (A.15) 

po 

Finally we treat region A. Here the correction reads: \/ zu'f -|- (1 — z)u2 — z{l — z) — 
zui — {l—z)u2, with the phase-space being Q({u+—2)'^ —u'^)@{—l+ zui+{l— z)u2— z{l— z)) , 
and 2 < u_|_ < oo and < u_ < 1. The numerical result is: 

r;^. = J-(_o.610Ca + 0.154Ca - 0.035n/). (A.16) 

PO 

Thus the final result for the non-inclusive correction factor r^j, which is + r^j -|- 
r^j + '^rai + '"ni multiplied by a factor 2, accounting for the fact that we have just considered 
U- > 0, is: 

2 

r„i = — (-2.145CA + O.GIOCa - 0.097n/) 

PO 

= -^(-3.071Ca - 0.194n/), (A.17) 

PO 

where we have written the result in all the above as a sum of three terms to show the 
contributions from the soft gluon, hard gluon and hard quark matrix elements respectively 
as in H^. 

Finally we note that f5o = 11/3Ca - 2/3n/ = 9(11) for nj = 3(0). 
A. 4 Comparison between kt and anti-kt results 

We present here a comparison between the r„j results for the kt and anti-^t results in 
various regions of the phase-space. A factor 2/(3q is left out: 
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Region 


anti-fcf 




A 

I 

B 
C 
D 


n ^'^7r^ . _i_ n 1 oq/^ . n ni fi>n ^ 

— U.00(L/y4 + U.iZoL/y4 — U.UiD?1.j 

-O.OOOC^ + O.OOOC^ - O.OOOn/ 
-0.245Ca + 0.0390 A - 0.004n/ 
-0.982Ca + O.SIOCa - 0.017n/ 
+0.557CA - 0.107CA - 0.015n/ 


— i.iDlL/yl + yj.Li iua — u.uoinj 
-O.OOSCa + O.OOICa - O.OOOn f 
-0.288Ca + 0.056Ca - O.OOGnj 
-1.2UCa + 0.382Ca - 0.025n/ 
+0.557Ca - O.mCA - 0.015n/ 


Sum 


-1.227Ca + 0.365Ca - 0.052n/ 


-2.U5Ca + O.GIOCa - 0.097n/ 



Table 3: Comparison between the non-inclusive coefficient r„i in the kt and anti-fct algorithms. 
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